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>P ■ Abstract. We study the consequences of superconducting quark cores (with color-flavor-locked phase as represen- 

tative example) for evolution of temperature profiles and the cooling curves in quark-hadron hybrid stars and in 
-..z , hypothetical self-bounded objects having no a hadron shell (quark core neutron stars). The quark gaps are varied 

<i ■ from to Ag = 50 MeV. For hybrid stars we find time scales of 1 -^ 5, 5 -^ 10 and 50 -ir 100 years for the formation of 

Oh' a quasistationary temperature distribution in the cases Ag = 0, 0.1 MeV and Z 1 MeV, respectively. These time 

scales are governed by the heat transport within quark cores for large diquark gaps (A iJ 1 MeV) and within the 
5_j ' hadron shell for small diquark gaps (A ;S 0.1 MeV). For quark core neutron stars we find a time scale ~ 300 years 

tzJ I for the formation of a quasistationary temperature distribution in the case A ;J 10 MeV and a very short one 

j^ . for A ;$ 1 MeV. If hot young compact objects will be observed they can be interpreted as manifestation of large 

gap color superconductivity. Depending on the size of the pairing gaps, the compact star takes different paths in 
the Ig(Ts) vs. lg(t) diagram where Ts is the surface temperature. Compared to the corresponding hadronic model 
^^ ■ which well fits existing data the model for the hybrid neutron star (with a large diquark gap) shows too fast 

Vh ' cooling. The same conclusion can be drawn for the corresponding self-bound objects. 

C^ \ 
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1. Introduction of one color whereas in the color-flavor locking (CFL) 

phase all the quarks are paired. 
The interiors of compact stars are considered as systems 
where high-density phases of strongly interacting matter Esti mates of the cooling evolution have been per 



do occur in nature, see Glendenning (1996) and Weber formed ( |Blaschke et al. 2000| ) for a self-bound isothermal 



(1999) for recent textbooks. The consequences of differ- quark core neutron star (QCNS) which has a crust but no 

ent phase transition scenarios for the cooling behaviour of hadron shell, and for a quark star (QS) which has neither 

compact stars have been reviewed recently in comparison crust nor hadron shell. It has been shown there in the case 

with existing X-ray data, see Page (1992), Schaab et al. of the 2SC (3SC) phase of QCNS that the consequences 

(1997). of the occurrence of gaps for the cooling curves are simi- 

A completely new situation might arise if the 1^^ *« the case of usual hadronic neutron stars (enhanced 

s cenarios suggest e d for (color) superconductivity cooling). However, for the CFL case it has been shown 

( lAlford et al. 19981 |Rapp et al. 199^ ) with large di- ^^^^ ^he cooling is extremely fast since the drop in the 



quark pairing gaps (A, - 50 -^ 100 MeV) in quark matter specific heat of superconducting quark matter dominates 

are applicable to neutron star interiors. Various phases o^cr the reduction of the neutrino emissivity. As has been 

are possible. The two-fiavor (2SC) or the three-flavor Pointed out there, the abnormal rate of the temperature 

(3SC) superconducting phases allow for unpaired quarks drop is the consequence of the approximation of homo- 

geneous temperature proflles the applicability of which 

Send offprint requests to: D. Blaschke should be limited by the heat transport effects. Page et 
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al. (2000) estimated the cooling of hybrid neutron stars 
(HNS) where heat transport effects within the supercon- 
ducting quark core have been disregarded. Neutrino mean 
free path in color superconducting quark matter have been 
discussed in ( Carter fc Reddy 2000| ) where a short period 
of cooling delay at the onset of color superconductivity 
for a QS has been conjectured in accordance with the es- 



timates of (Blaschke et al. 2000) in the CFL case for small 
gaps. 

In the present paper we want to consider a more de- 
tailed scenario for the thermal evolution of HNS and 
QCNS which includes the heat transport in both the quark 
and the hadronic matter. We will demonstrate how long it 
takes for the HNS and for the QCNS to establish a quasis- 
tationary temperature profile and then we consider the in- 
fluence of both the diquark pairing gaps and the hadronic 
gaps on the evolution of the surface temperature. 

2. Processes in HNS and QCNS 

2.1. Processes in quark matter 
2.1.1. Emissivity 

A detailed discussion of the neutrino emissivity of quark 
matter without the possibility of color superconductivity 
has been given first by Iwamoto (1982). In his work the 
quark direct Urea (QDU) reactions d -^ uev and ue -^ dv 
have been suggested as the most efficient processes. Their 
emissivities have been obtained as n 

e?°" ~ 9.4 X 1026a,ur;/3^QDu T^ erg cm'^ ^-i^ ^^^ 

where at a compression u = n/riQ ~ 2 the strong cou- 
pling constant is a^ « 1 and decreases logarithmically 
at still higher densities. The nuclear saturation density 
is no = 0.17 fm~ , Yg = rie/n is the electron fraction, 
and Tg is the temperature in units of 10^ K. If for some- 
what larger density the electron fraction was too small 
(Ye < Y(ic — 10~*), then a ll the QDU process es would 
be completely switched off ( Duncan et al. 1983 ) and the 
neutrino emission would be governed by two-quark re- 
actions like the quark modified Urea (QMU) and the 
quark bremsstrahlung (QB) processes dq —f uqeD and 
9192 ~> 9192^1^, respectively. The emissivities of the QMU 
and QB processes have been estimated as ([wamoto 1982) 



^QMU 



e^^ ~ 9.0 X IO^^Cqmu T| erg cm-^ g-i. 



(2) 



Due to the pairing, the emissivities of QDU processes 
are suppressed by a factor Cqdu ~ exp(— A^/T) and 
the emissivities of QMU and QB processes are sup- 
pressed by a factor Cqmu ~ exp(— 2Aq/T) for T < 
T'crit.g — 0.4 Aq whcrcas for T > Tcrit.g these factors 
are equal to unity. The modification of Tcrit,q{Aq) rel- 
ative to the standard BCS formula is due to the for- 
mation of correlations as, e.g., instanton- anti-instanton 

^ Please notice that the numerical factors in all the esti- 
mates below depend essentially on rather unknown factors and 
thereby can be still varied. 



molecules ( Rapp et al. 200C| ). For the temperature depen- 
dence of the gap below Td-n^q we use the interpolation 
formula A(r) = A{Q)y/l -T/T„it,q, with A(0) being the 
gap at zero temperature. 

The contribution of the reaction ee — > eevD is very 
small ( Kaminker fc Haensel 1999| ) 

(3) 



eff = 2.8 X 10'^ r 



1/3,1/3^8 



U 



-3 „-l 



Tg erg cm ■^ s 



but it can become important when quark processes are 
blocked out for large values of Aq/T in superconducting 
quark matter. 

2.1.2. Specific heat 

For the quark specific heat we use an expression 
([Iwamoto 1982D 



Cq ~ 10^1^2/3^3 Tg erg cm'^ K~\ 



(4) 

where Cs - 3.1 (Terit,,/^)^/^ exp(-A,/r). Besides, 
one should add the gluon-photon contribution 



(5) 



( Plaschke et al. 2000D 

Cg-^ = 3.0 X 10^^ Ng^j T| erg cm'^ K'^, 

where A^g-7 is the number of available massless gluon- 
photon states (which are present even in the color super- 
conducting phase), as well as the electron one, 

Ce = 5.7 X 10^^ Y^/^u^/^ Tg erg cm"^ K'\ (6) 



2.1.3. Heat conductivity 

The heat cond uctivity of the matte r is the sum of partial 
contributions ( Flowers fc Itoh 1981 ) 

' -^^ P) 



, l^ii 



i 



where i,j denote the components (particle species). 
For quark matter k is the sum of the partial conduc- 
tivities of the electron, quark and gluon components 



(8) 



( [Haensel fc Jerzak 1989| , [Baiko fc Haensel 1999D 

where k^ ~ K^e is determined by electron-electron scatter- 
ing processes since in superconducting quark matter the 
partial contribution 1/ n^q (as well as ^/ K-gq ) is addition- 
ally suppressed by a Cqdu factor, as for the scattering on 
impurities in metallic superconductors. For Kee we have 

Kee = 5.5 X lO^^u Ye Tg~^ erg s^^cm^^K^^ 

and 



^qq 



(9) 



~ 1.1 X 10^3 J— u CQDuTg-^erg s-^cm-^K-i, (10) 

V "s 

where we have accounted for the suppression factor. The 
contribution of massless gluons we estimate as 



^g — '^gg 

~ 6.0 X 10" T^ erg s-^cm^^K-^ 



(11) 
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2.2. Processes in hadronic matter 



2.2.1. Emissivity 



gap which we took as function of density from Fig. 2 of 
( ISchaab et al. 1997| ). 

For T < Tcrit the most important contribu- 
DU processes n -^ peD, pe -^ nv can occur when the pro- tion comes from the neutr on ([Flowers et al. 1976 , 



ton fraction exceeds 11% (Lattimcr et al. 1991). It does 
not occur for the equation of state wc use but there are 



examples for the opposite case (Page et al. 2000). This is 
another difference of our present work to the recent one 
of Page et al. (2000). 

Next, we take into account the modified Urea 
(MU) processes nn — > npev, np -^ ppev, and the 
r everse processes. The final emissivity is given by 
(iFriman fc Maxwell 1979| , [Yakovlev et al. 1999|) 



j;^u = 8.6 X 102im;tiu(reu)i/3c„MuT«ergcm-3s-i(12) 
P^u ^ 8.5 X 102im;i,u(>;")'/'CpMur«ergcm-3s-i(13) 



Here m* = \/m*^^ i + Pp i ^^ ^^^ non-relativistic quasi- 
particle effective mass related to the in-medium one- 
particle energies from a given relativistic mean field 
model for i — n,p. We have introduced the abbre- 



viations m 



nMU 



= {m*Jmn) {m*/mp) and m. 



■pMU 



(nip/mp) (m'^/mn). The suppression factors are CnMU ~ 
CnCp ^ exp{-[A„(r)+Ap(T)]/T}, CpMU ^ Q, and should 
be replaced by unity for T > Tcrit, i when for given species 
i the corresponding gap vanishes. For neutron and proton 
5'-wave pairing, Ai(0) = 1.76 Tcrit, i and for P-wave pair- 
ing of neutrons A„(0) = 1.19 Tcrit, «. To be conservative 
we have used in (02) the free one-pion exchange estimate 
of the NN interaction amplitude. Restricting ourselves 
to a qualitative analysis we use here simplified exponen- 
tial suppression factors C,i. In a more detailed analysis 
these Ci-factors have prefactors with rathe r strong tem- 
perature dependences Yakovlev et al. 1999 . At tempera- 
tures T ^ Tc their inclusion only slightly affects the re- 
sulting cooling curves. For T <^ Tc the MU process gives 
in any case a negligible contribution to the total emissiv- 
ity and thereby corresponding modifications can again be 
omitted. Also for the sake of simplicity the general pos- 
sibility of a ■^P2(|m^j| — 2) pairing which may result in 
a power-law behaviour of the specific heat and the emis- 



sivity of the MU process Voskresensky & Senatorov 1987, 
Yakovlev et al. 1999 is disregarded since mechanisms of 
this type of pairing are up to now not elaborated. Even 
more essential modifications of the MU rates may pre- 
sumably come from in-medium effects which could result 
in extra prefactors of 10^ -^ 10^ already at T ~ T^. 

In order to estimate the role of the in-medium effects 
in the NN interaction for the HNS cooling we have also 
performed calculations for the so-called medium modified 
Urea (MMU) process (Voskresensky & Senatorov 1986, 



Migdal et al. 1990| ) by multiplying the rates (|lj) by the 
appropriate prefactor 



^MMU /^MU 



/e^''^l0^r'{g')/^'{kc,pF)]u''/^ 



(14) 



where the value r{g') ~ 1/[1 -I- 1.4 m^/^] is due to the 
dressing of irNN vertices and uj < m-^ is the effective pion 



Voskresensky fc Senatorov 198'i , Yakovlev et al. 1999| ) 



and the proton ( Voskresensky fc Senatorov 1987 ) pair 
breaking and formation processes. Their emissivities we 
take from Ref. (Schaab et al. 1997) which is applicable 
for both the cases of S- and P-wave nucleon pairing 



,NPBF 



,PPBF 



where 



6.6 X 102^(m;/TO„)(A„(T)/MeV)^ u^/^ 

xe/(A„(T)/T)ergcm-V\ 

0.8 X 102*^(m*/mp)(Ap(T)/MeV)^ w^/^ 



-3^-1 



X /(Ap(T)/T) erg cm^-'s 



(15) 



(16) 



/(A,(T)/T) ~ 0.89 Vr/A,(T)exp[-2A,(T)/T] , (17) 

f ~ 0.5 for ^5o pairing and ^ ~ 1 for ^Pj pair- 
ing. A significant contribution of the proton channel is 
due to the NN correlation effect taken into account in 



(Voskresensky fc Senatorov 1987) 



2.2.2. Specific heat 

For the nucleons {i 
( [Maxwell 1979| ) 



n,p), the specific heat is 



Q = 1.6 X 10'*'(m*/m,) u^/^^mmu Tg erg cm^^j^-^ (18) 

The specific heat for the electrons is determined by (pi) 
and that for the photons for T > Tcp is given by (H) with 
the corresponding number of polarizations A'-y — 2. 



2.2.3. Heat conductivity 

The total conductivity is 



(19) 



where the partial contributions are (Flowers fc Itoh 1981, 
Baiko fc Haensel 1999D 



^nn 


= 8.3 X 10^2 f™"' 


4 S A 
Skn Tq 


Skn 


= 0.38 z^^/2 + 3.7 


,2/5 


KjfLp 


= 8.9x10^6 (""""' 


2 ,2/- rp 

1 Zpkp 



erg s"^cm-iK"\ (20) 

(21) 

%^ erg s-^cm-iK-i, (22) 

Skp = 1.83 z^2 + 1.43 ^2(Q 4 _^ ^8)^1 (23) 

Here we have introduced the apropriate suppression fac- 
tors Ci which act in the presence of gaps for superfluid 
hadronic matter and we have used the abbreviation Zi = 
(ni/(4no))^/'^. The heat conductivity of electrons is given 
by Eq. (I). 
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2.3. Processes in the mixed piiase 

In some density interval nn < n < uq we have a mixed 
phase in the HNS case. Although our treatment of the 
contribution of the mixed phase to the equation of state 
is rather poor, for the discussion of the neutrino pro- 
cesses we assume the mixed phase to be constructed as 
a lattice (for T < Tmoit) of color superconducting quark 
droplets embedded in the nucleonic Fermi sea (superfluid 
for T < Tcrit.i; i = n,p) at densities near the lower bound- 
ary density rifj which is reorganized as a lattice of nucleon 
droplets in the quark Fermi sea near the upper boundary 
density rtq . We suppose that the mixed phase contributes 
to the neutrino reactions and to the heat transport as two 
homogeneous phases weighted according to their partial 
volumina and thus suppress reactions on the impurities 
which could be as efficient as one-nucleon/quark processes 
in the liquid phase (Reddy et al. 200C). Thus we assume 
T < Tincit in our case. 



2.4. Crust and surface 

In order to avoid a complicated analysis of the pro- 
cesses in the crust we use an interpolated relation 
between the surface and the crust temperatures 
(see Shapiro fc Tcukolsky 198^ ) T, = (10 T,„)2/3, 
where temperatures are given in K, the so- 
called Tsuruta law Tsuruta 1979. More elabo- 



rated and complicated dependences of Tg on T™ 
and the parameters of the crust can be found in 



(Gudmundsson et al. 1982, Gudmundsson et al. 1982 



Van Riper 198? , Potekhin et al. 1997). As we have es- 
timated, their inclusion does not affect qualitatively 
the conclusions of this work. The mantle temperature 
T„i — T{rm) we take at the distance r^ from the 
center where the density is equal to the critical one 
for the transition to the mixed (Aen) phase. This is 
correct for time scales exceeding those determined by 
the heat conductivity of the crust estimated as ^ 10 
yr for hadronic stars and ^ 1 yr for strange stars (see 
Pizzochero 1991 , Lattimcr et al. 1991 ). Since the crust is 



rather thin and has a much smaller heat content than the 
bulk of the star its contribution to the cooling delay for 
time intervals shorter than the ones given above is negli- 
gible. Due to these reasons and for the sake of simplicity 
we disregard these effects in further considerations. 

We also add the photon contribution to the cmissivity 
at the surface 



2 X 10 



18 



R 



10 km 



Tg-j erg cm ^ sec ^, 



(24) 



where Tgj is the surface temperature in units of 10^ K, 
using it as the boundary condition in our transport code. 
Thus simplifying we assume that the outer crust and the 
photosphere are thin enough to approximate r^ ~ i?, 
where R is the star radius. Then, the luminosity at r,„ 
corresponds to the surface luminosity of the photons from 
the star /(r^) = L^{Ts). 



3. Evolution of inhomogeneous emissivity profiles 

3.1. Equation of state and gaps 

The equation of state is decribed according to the 



model (Chubarian et al. 2000) incorporating three possi 



ble phases (depending on the star mass, the value of the 
bag constant B and other parameters as the mass of the 
strange quark, the compressibility of hadron matter etc.). 
These phases are hadronic matter, quark-hadron mixed 
phase and pure quark matter. As a representative exam- 
ple we consider below a star with a mass of 1.4 Mq for 
the neutron star consisting of pure hadronic matter (ra- 
dius ~ 13 km), for the HNS (radius ~ 10 km) and for the 
self-bound QCNS (radius ~ 9 km). The typical electron 
fraction is Fg ~ 0.04 for the hadronic star which does not 
allow for the occurrence of direct Urea processes. 

For the HNS the electron fraction is Yg ~ 0.03 for 
n = TiH , where tih corresponds to the critical density when 
first seeds of quark phase appear; Y^ ~ 8.5 x 10~^ for 
n = TiQ when hadron seeds disappear, i.e. at the boundary 
of the quark core, and Yg — 2 x 10^^ for n = ric, i.e. for 
the baryon density in the star center. 

The radius corresponding to n = nu is ru = 8 km and 
that corresponding to n = nq is rq = 6.5 km. We take 
the same dependences Ai(n) for i = n,p as in Fig. 6 of 
(Schaab et al. 1997). As an alternative example we also 



consider the model with suppressed gaps to estimate the 
effect of possible uncertainties. In the mixed phase we take 
Ai = Ai(n-ti) for the hadron sub-system and A = A, for 
the quark one, and in the quark phase we for simplicity 
take Aq = const. 

For QCNS besides the case when Y^ is determined from 
our calculation we consider also Yg = as an extreme 
example of the case Yg < Y'gg, in order to understand its 
possible consequencies. Then direct Urea does not occur 
also in the quark phase. 



3.2. Cooling 

Due to the inhomogeneous distribution of matter inside 
the star and finite heat conductivity the temperature pro- 
file T(r, t) during the cooling process can differ from the 
isothermal equilibrium one for which the temperature on 
the inner crust boundary (T^) and the central one (Tg) 
are connected by the relation T„i — Tgexp[<I>(0) — $(i?)], 
where $(0) — ^{R) is the difference of the gravitational 
potentials in the center and at the surface of the star, 
respectively. 

The flux of energy l{r) per unit time through a spher- 
ical slice at the distance r from the center is proportional 
to the gradient of the temperature on both sides of the 
spherical slice. 



l{r) — —Attt K{r) 



9(re*) 
dr 



1 



2M 



(25) 



the factor e * W 1 — 2M corresponds to the relativistic 
correction of the time scale and the unit of thickness. The 
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equations for energy balance and thermal energy trans- 
port are ( [Weber 1999| ) 



M= 1.4M 



_d_ 
dA 



(/e-) 



— fTe^*) 



ot 



(26) 



(27) 



where n = n{r) is the baryon number density, A = A(r) is 
the total baryon number in the sphere with radius r and 



dr 

dA 



1 



A'Kr'^Ti 



1 



2Af 



The total neutrino emissivity e^, and the total specific heat 
cv are given as the sum of the corresponding partial con- 
tributions defined in the previous Section for a composi- 
tion riiir) of constituents i of the matter under the condi- 
tions of the actual temperature profile T{r,t). The accu- 
mulated mass M — M(r) and the gravitational potential 
$ = <i>(r) can be determined by 



dM 
'dA 

dA 



1 



Anr^p - 



2M 

r 
- m 



Airr'^j 



1 



2 A/ 



(29) 
(30) 



where e = e{r) is the energy density profile and the pres- 
sure profile p — p{r) is defined by the condition of hydro- 
dynamical equilibrium 



dp , ,d^ 

5^ = -(P + s)^ 



(31) 



The boundary conditions for the solution of ( pq ) and (E^ 
read l{r = 0) = /(A = 0) = and T(A(r„) = A,t) = 
T{rm ^ R,t) ^ T„i{t), respectively. 

In our examples we choose the initial temperature to 
be 1 MeV. This is a typical value for the temperature 
2^opacity st which the star becomes transparent for neu- 
trinos. Simplifying we disregard the neutrino infiuence on 
transport. These effects dominate for i < 1 -^ 100 min, 
when the star cools down to T < Topacity and become 
unimportant for later times. 

4. Temperature evolution of HNS 

4.1. Temperature profiles T{r, t) 

With the above inputs we have solved the evolution equa- 
tion for the temperature profile. In order to demonstrate 
the influence of the size of the diquark and nucleon pairing 
gaps on the evolution of the temperature profile we have 
performed solutions with different values of the quark and 
nucleon gaps. The representative examples are shown in 
Fig. 1. Nucleon gap s are taken the same as in Fig. 6 of 
( ^chaab et al. 1997 ). The quark gaps larger than 1 MeV 
show the same typical behaviour as for the gap Ag — 50 
MeV. The gaps much smaller than 1 MeV exhibit the 
typical small gap behaviour as for A^ = 0.1 McV in our 
example. 



o 



A^:\&X 



^ 0.1 

s 

^ 0.01 



(28) S 



> 

I 0.1 



0.01 



0.01 



.'y 



I I .- 



t = 
10"^ yr 
0.1 yr 
lyr 
5 yr 
50 yr 
100 yr 
500 yr 



A =0.1 MeV 

q 




. _ A = 50 MeV 



_L 



0.5 



r/r 



Fig. 1. Early evolution of temperature profiles for a HNS 
of M = \.\Mq with large (lower panel), small (middle) 
and vanishing (upper panel) diquark pairing gaps. 



For large quark gaps, the hadron phase cools down on 
the time scale i ~ 0.1 -^ 5 yr due to the heat transport to 
the star surface whereas the quark core keeps the heat dur- 
ing this time. Therefore, if an independent measurement 
of the core temperature was possible, e.g. by neutrinos, 
then a core temperature stall during first several years of 
cooling evolution would be a case for quark core super- 
conductivity with large pairing gaps. Then the quark core 
begins to cool down slowly from its mixed phase bound- 
ary, the cold from hadrons via mixed phase spreads to the 
center, demonstrating that the direct neutrino radiations 
during all the time are unefficient within the quark core. 
The homogeneous temperature profile is recovered at typ- 
ical times of i ~ 50 -^ 100 yr. 

For Aq = 0.1 MeV the QDU processes and the heat 
conductivity within the quark core are quite efficient. 
Therefore at very small times (i <C 10~^ yr) the tem- 
perature profile within the quark core is a homogeneous 
one whereas during the next 5 -^ 10 yr the hadron shell is 
cooled down from the both sides, the mixed phase bound- 
ary at n(rH) = nn and the star surface. Thus, for A^ ^ 1 
MeV the interior evolution of the temperature has an in- 
fluence on the surface temperature from the very early 
times. Actually the shortest time scale which determines 
the slowing of the heat transport within the crust is in 
this case i ~ 1 yr. Thus at much smaller times one has 
Ts{t) ~ Ts(0) with a drop of T(r, i) in the vicinity of 
T JTm — 1 . Due to the relatively small heat content of the 
crust this peculiar effect is disregarded here. 
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Qualitatively the same behaviour is illustrated for nor- 
mal quark matter with the only difference that quantita- 
tively the quark processes are more efficient at A^ = and 
the cold thereby traverses the hadron phase more rapidly, 
at typical times i ~ 1 -^ 5 yr. 



4.2. Evolution of the surface temperature 

A detailed comparison of the cooling evolution (Ig Ts vs. 
Ig t) of HNS for different values of quark and hadron 
gaps is given in Fig. 2. We have found that the curves 
for Ag > 1 MeV are very close to each other demonstrat- 
ing typical large gap behaviour. As representative example 
we again take Ag = 50 MeV. The behaviour of the cool- 
ing curve for i < 50 -^ 100 yr is in straight correspondence 
with the heat transport processes discussed above. The 
subsequent time evolution is governed by the processes in 
the hadronic shell and by a delayed transport within the 
quark core with a dramatically suppressed neutrino emis- 
sivity from the color superconducting region. In order to 
demonstrate this feature we have performed a calculation 
with the nucleon gaps (Ai(n), i = n,p) being artificially 
suppressed by a factor 0.2. Then up to lg(i[yr]) < 4 the 
behaviour of the cooling curve is analogous to the one we 
would have obtained for pure hadronic matter. The curves 
labelled " MMU" are calculated with the rates of modified 
Urea processes of Eq. ( [l2| ) multiphed by the factor (|l4|), 
i.e. with inclusion of appropriate medium modifications in 
the NN interaction. As can be seen from Fig. 2, these ef- 
fects have an infiuence on the cooling evolution only for 
lg(f [yr]) < 2 since our specific model equation of state does 
not allow for large nucleon densities in the hadron phase 
for the M = 1.4 Mq neutron star under discussion. The 
effect would be more pronounced for larger star masses, 
a softer equaton of state for hadron matter and smaller 
values of the gaps in the hadronic phase. 

The unique asymptotic behaviour at lg(i[yr]) > 5 for 
all the curves corresponding to finite values of the quark 
and nucleon gaps is due to a competition between normal 
electron contribution to the specific heat and the photon 
emissivity from the surface since small exponentials switch 
off all the processes related to paired particles. This tail is 
very sensitive to the interpolation law T^ — f{Tm) used to 
simplify the consideration of the crust. The curves coincide 
at large times due to the uniquely chosen relation T^ oc 

rp2/3. 

The curves for A^ = 0.1 MeV demonstrate an interme- 
diate cooling behaviour between those for Ag = 50 MeV 
and Ag = 0. Heat transport becomes not efficient after 
first 5-^10 yr. The subsequent 10^ yr evolution is governed 
by QDU processes and quark specific heat being only mod- 
erately suppressed by the gaps and by the rates of NPBF 
processes in the hadronic matter (up to lg(t[yr]) < 2.5). 
At lg(i[yr]) > 4 begins the photon cooling era. 

The curves for normal quark matter (A^ = 0) are gov- 
erned by the heat transport at times i ;$ 5 yr and then by 
QDU processes and the quark specific heat. The NPBF 
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log(t[yr]) 
Fig. 2. Evolution of the surface temperature Ts of HNS 

2/3 

with M ~ 1.4i\/Q for T^ ex Tm . Data points are from 
Schaab et al. (1999) (full symbols) and from Yakovlev et 
al. (2000) (empty symbols), tw is the typical time which is 
necessary for the cooling wave to pass through the crust. 

processes are important up to lg(i[yr]) < 2, the photon 
era is delayed up to lg(i[yr]) > 7. For times smaller than 
tw (see Fig. 2) the heat transport is delayed within the 
crust area (Lattimer et al. 1991). Since we, for simplicity, 
disregarded this delay in our heat transport treatment, for 
such small times the curves should be interpreted as the 
Tm{t) dependence scaled to guide the eye by the same law 
oc J m , as J o . 



5. Temperature evolution of QCNS 

5.1. Temperature profiles T{r, t) 

As can be seen from the lower panel of Fig. 3 for our 
large gap example (A, = 50MeV), where l^e 7^ is de- 
termined from equation of state, the QCNS cools down 
rather smoothly from its surface during the first 300 yr. 

All the small gap examples exhibit about the same be- 
haviour as the case A^ = 0.1 MeV which we present in the 
upper panel of Fig. 3. Opposite to the large gap case, the 
transport processes are very efficient and a homogeneous 
temperature profile is recovered at t ^ 10~^ yr. In real- 
ity in the latter case T{rm,t) ~ T(rm,0) at times ^ 1 yr 
with a drop of T(r) in a narrow interval near r/rm — 1 and 
we disregard this peculiar behaviour as we have explained 
this above discussing HNS. 
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QCNS with M = 1.4A/q for T^ oc T, 



5.2. Evolution of the surface temperature 

In Fig. 4 we demonstrate the evolution of the surface tem- 
perature Ig Ts vs. Ig t for representative cases. The curve 
for Aq = 50 MeV {Yf, ^ 0) shows a 300 yr delayed cooling 
which then is controlled by the electron specific heat and 



the photon emissivity from the surface, as it was stated 
in ( Blaschke ct al. 200C| ). However, in difference with the 
corresponding curve of (Blaschke et al. 200C) we have a 
much larger typical cooling time scale (i ~ 300 yr) since 
in that work the simplifying approximation of a homoge- 
neous temperature profile has been used. For the low gap 
limit Ag ^ 10 MeV in this particular QCNS case there 
is no heat transport delay and our results coincide with 



those obtained in (Blaschke et al. 2000). Unimportant dif- 
ferences are only due to the inhomogeneous density profile 
and a different value of Ye{n) which follows in our case 
from an actual calculation for the given equation of state. 
The curve corresponding to Ag = 0.1 MeV (Yg 7^ 0) from 
very early times is due to a competition between QDU 
emissivity and the quark specific heat, both being only 
moderately suppressed whereas late time asymptotics, as 
in the large gap example, relates to the photon era. A 
similar trend is present for normal quark matter staying 
below the corresponding curves for Ag = 0.1 MeV due 
to the absence of the C suppression factors in the former 
case. At t > 10^ yr both cases have a different photon era 
asymptotic behaviour whereas all finite gap curves have 
the same large time asymptotics when Ye ^ 0. 

At Yg = the QDU processes are absent and the large 
time asymptotics are governed by the giuon-photon spe- 
cific heat which is a nonlinear function of the tempera- 
ture and the photon emissivity. Due to this nonlinearity 
the large time asymptotics are different in all the cases 
under consideration. All the results for Ye — (includ- 
ing case of large gaps) coincide with those obtained in 
(Blaschke et al. 2000) although the sharp fall of the curve 
for the large gap case is in reality controlled by the heat 
conductivity of the crust that will lead to a cooling delay 
within the first year of evolution ( Pizzochero 1991 ) as in 
Fig. 2 the curves for t <tw should be interpreted as Tm(t) 
dependence scaled by the Tsuruta law. 



6. Conclusion 

For the CFL phase with large quark gap, which was ex- 
pected to exhibit the most prominent manifestations of 
colour superconductivity in HNS and QCNS, we have 
found an essential delay of the cooling during the first 
50 ^ 300 yr due to a dramatic suppression of the heat con- 
ductivity in the quark matter region. This delay makes the 
cooling of QCNS not as rapid as one could expect when 
ignoring the heat transport. In HNS compared to QCNS 
(large gaps) an additional delay of the subsequent cool- 
ing evolution comes from the processes in pure hadronic 
matter. In spite of that we found still too fast cooling 
for those objects. Therefore, with the CFL phase of large 
quark gap it seems rather difficult to explain the major- 
ity of the presently known data both in the cases of the 
HNS and QCNS, whereas in the case of pure hadronic 
stars the available data are much better fitted within the 
same model for the hadronic matter we used here. For 
2SC (3SC) phases we expect analogous behaviour to that 
demonstrated by Ag = since QDU processes on unpaired 
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quarks are then allowed. We however do not exclude that 
new observations may lead to lower surface temperatures 
for some supernova remnants and will be better consistent 
with the model which also needs further improvements. 
On the other hand, if future observations will show very 
large temperatures for young compact stars they could be 
interpreted as a manifestation of large gap color supercon- 
ductivity in the interiors of these objects. 
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